Detection and control of individual nuclear spins using a weakly coupled electron spin 
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We experimentally isolate, characterize and coherently control up to six individual nuclear spins 
that are weakly coupled to an electron spin in diamond. Our method employs multi-pulse sequences 
on the electron spin that resonantly amplify the interaction with a selected nuclear spin and at the 
same time dynamically suppress decoherence caused by the rest of the spin bath. We are able to 
address nuclear spins with interaction strengths that are an order of magnitude smaller than the 
electron spin dephasing rate. Our results provide a route towards tomography with single-nuclear- 
spin sensitivity and greatly extend the number of available quantum bits for quantum information 
processing in diamond. 



Detecting the weak magnetic moment of a single nu- 
clear spin presents the ultimate Umit of sensitivity in 
magnetic resonance imaging [TH3]. Furthermore, nuclear 
spins may play a key role as qubits with long coherence 
times in quantum information technologies Address- 
ing and controlling single nuclear spins is challenging be- 
cause the spins are generally embedded in a noisy envi- 
ronment, such as a surrounding bath of nuclear spins. 

The electron spin of a nitrogen-vacancy (NV) cen- 
ter is a powerful probe of its local magnetic environ- 
ment [H [21 [SHU] • If a single or a few nuclear spins are 
located particularly close to an NV center, the hyper- 
fine interaction can well exceed the electron spin dephas- 
ing rate I/T2* pjj. Such strongly coupled nuclear spins 
are readily distinguished from the rest of the spin bath 
[Ol [13] and can be selectively addressed and controlled 
[T51 - [2^ . However, typically the nuclear spin of interest 
is embedded in a bath of fluctuating nuclear spins. As 
a result, the coupling of this single nuclear spin to the 
NV center is weak compared to the rate of electron spin 
dephasing induced by the spin bath. For both magne- 
tometry and quantum information purposes it would be 
greatly beneficial to be able to individually resolve and 
address such weakly coupled nuclear spins. 

In this Letter, we isolate, characterize and selectively 
control up to six weakly coupled ^^C nuclear spins that 
are embedded in the spin bath surrounding an NV cen- 
ter. The weak signal of a specific nuclear spin is amplified 
by precisely tuning a multi-pulse control sequence on the 
NV electron spin into resonance with the electron-nuclear 
spin dynamics [23]. At the same time this sequence dy- 
namically decouples the electron spin from all other nu- 
clear spins [23H2S]- With this technique, we are able to 
resolve and coherently control nuclear spins with cou- 
plings that are an order of magnitude smaller than the 
dephasing rate of the NV center. Our results can enable 
tomography with single nuclear spin sensitivity and have 
the potential to greatly extend the number of solid-state 



spin qubits available for quantum information processing. 

Our method to isolate a weakly coupled nuclear spin 
from a background of other nuclear spins is based on 
the distinct conditional precession of each nuclear spin 
due to its particular hyperfine interaction with the NV 
electron spin {S = 1), Fig. 1(a). For the electron in nis = 
0, all nuclear spins precess with the Larmor frequency 



(71 71 \l 

r I 2r I r ) 



N/2 




m = m = -1 



m = m = -1 



FIG. 1; Concept of isolating and controlling weakly cou- 
pled spins, (a) Surrounding ^^C nuclear spins precess about 
axes that depend on the NV electron spin state. For rUs — 0, 
all ^^C spins precess about ljl set by the applied magnetic 
field Bq. For rua = — 1, each spin precesses about a distinct 
axes a) due to the hyperfine interaction u^h- (b) Calculated 
probability to preserve the initial electron spin state af- 
ter a decoupling sequence with A'^ = 32, for two ^'^C spins 
with e = 7r/4.5, Nucleus I: a;^ = 27r ■ 40 kHz, Nucleus II: 
ojh = 27r • 20 kHz, _Bo=293 C. Each spin can be selectively 
addressed by tuning the inter-pulse delay 2r into resonance 
with its dynamics. (c,d) Bloch spheres showing the nuclear 
spin dynamics for r resonant with nucleus I (arrow), (c) For 
nucleus I, the net result is a rotation around anti-parallel axes 
(rio and rfi) for the two electron states, resulting in entangle- 
ment, (d) Nucleus II is decoupled: its rotation is independent 
of the electron state. 
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ijJl around an axis parallel to the applied magnetic field 
Bq. For vTLs = — 1, each nuclear spin precesses around a 
distinct axis Cj = ojl +^h- The hyperfine interaction a;^ 
depends on the position of that particular nuclear spin 
relative to the NV center. 

We can probe this conditional interaction by preparing 
the electron spin in a superposition, \x) — {\ms = 0) + 
\ms = — 1))/-\/2, and applying a dynamical decoupling 
sequence consisting of N sequential 7r-pulses. Consider 
the basic decoupling unit on the electron spin r — tt — 2r — 
TT — T, in which r is a free evolution time (Fig. 1(b)). The 
net result of this unit is a rotation of the nuclear spin by 
an angle (j) around an axis that depends on the initial 
state of the electron spin: no for initial state nis = and 
fii for initial state nis = —1 ^51 [?7] . 

If no and rii are not parallel, the resulting conditional 
rotation of the nuclear spin generally entangles the elec- 
tron and nuclear spins. As a result, for an unpolarized 
nuclear spin state, the final electron spin state is a statis- 
tical mixture of \x) and \—x) = (|ms=0) — |ms=— 1))/-\/2. 
The probability that the initial state |a;) is preserved is 
given by 

= (M + l)/2, (1) 
with for a single nuclear spin: 

A/ = l-(l-no-ni)sin2^. (2) 

For multiple nuclear spins that do not mutually interact, 
M is given by the product of all the individual values Mj 
for each individual spin j. Analytical expressions for 
and for the angle between no and fii as a function of the 
hyperfine interaction uJh and the inter-pulse delay r are 
given in the supplemental material |27) . 

As an example. Fig. 1(b) shows calculated results for 
two ^'^C spins with different hyperfine interactions. For 
most values of r the NV spin is effectively decoupled from 
both nuclear spins and its initial state is conserved {Px ~ 
1). For specific values of r, the sequence is precisely 
resonant for one of the ^"^C spins and a sharp dip in the 
signal is observed. Figures l(c, d) illustrate the evolution 
of the nuclear spins at the resonance condition for nuclear 
spin I. At this value of r, the net rotation axes fio and 
rii for nuclear spin I are approximately an^i-parallel and 
the resulting conditional rotation entangles nuclear spin 
I with the electron spin (P^ « 1/2). In contrast, at the 
same value of r, no and fii are nearly parallel for nuclear 
spin II and the resulting unconditional rotation leaves 
the electron spin unaffected. These resonances appear 
periodically as a function of r. 

More insight into the periodicity and depth of the res- 
onances can be gained by considering the case of large 
magnetic field, a;^ 3> w/j. In this case the positions of 
the resonances are given by [27]: 



where k = 1,2,3.. is the order of the resonance, and 
A is the parallel component of the hyperfine interaction 
A = ujh cos 6. Equation |3] shows that the position is a 
linear function of k. The amplitude of the resonances is 
governed by the rotation angle (j), which is of order B/lul, 
with B = ujh sin 9 the perpendicular component of the 
hyperfine coupling. Although (j) is small, the total angle 
is amplified by the large number of pulses N, enabling the 
detection with maximum contrast even of weakly coupled 
spins. In this way a single nuclear spin can be isolated 
from a bath of spins by a judicious choice of the inter- 
pulse delay 2r and the number of pulses N. 

We experimentally demonstrate our method using an 
NV center in a type Ila diamond with a natural abun- 
dance of ^■^C nuclear spins (1.1%). All experiments are 
performed at room temperature with an applied mag- 
netic field along the NV symmetry axis. The NV elec- 
tron spin is prepared in = by illumination with a 
532 nm laser and read out through its spin-dependent flu- 
orescence. The experimental setup is described in detail 
in Ref. 23. 

We choose an NV center that shows no nearby strongly 
coupled ^'^C spins in the electron spin resonance (ESR) 
spectrum and Ramsey measurements. The hyperfine 
coupling to the NV spin of all individual ^■^C spins is thus 
weak compared to I/T2*: all individual nuclear spins are 
hidden in the spin bath. 

The experimental signal for a decoupling sequence with 
32 TT-pulses is shown in Fig. 2(a). We observe sharp dips 
and broader collapses in an approximately exponentially 
decaying signal (see Fig. 2(b) for a magnification). The 
broader collapses correspond to the overlapping signals 
of multiple nuclear spins in the spin bath, whose product 
tends to yield Px « 0.5. The sharp dips are signatures 
of the resonances of individual ^"^C spins. These appear 
primarily for large r because the separation between res- 
onances of different spins increases with the resonance 
order k (see Eq. |3]). We exploit the linear dependence 
in Eq. [3] to identify five distinct ^■^C spins (Fig. 2(c)). 
The resonances assigned to these spins are indicated in 
Fig. 2(b). 

With a fit based on Eq. [2] we are able to determine 
both the magnitude uj^ and the angle 9 of the hyperfine 
coupling from the experimentally observed resonances in 
Fig. 2(b) for each of the five spins. These fits take the 
overall signal decay due to relaxation to rris = -1-1 and de- 
phasing of the electron state into account [57] . Although 
nuclear spin 6 can not be clearly resolved from the spin 
bath with a sequence of 32 pulses (Fig. 2(b)), we can 
further increase the sensitivity by applying more pulses. 
For = 96 the signal for spin 6 is well-isolated from 
the spin bath (Fig. 2(d)), enabling the characterization 
of the hyperfine interaction. 

The obtained values for the hyperfine interaction 
strength Uh and angle 9 for the six nuclear spins are listed 
in Table 1. These values should be compared to the min- 
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FIG. 2; Resolving individual weakly coupled '^^C nuclear spins, (a) Px as function of r for a decoupling sequence with A'^ = 32 
and a magnetic field Bq = 401 G. The sharp resonances in the echo signal correspond to the coherent interaction with individual 
^^C atoms, (b) Magnification of the section marked in (a) indicating resonances associated with six nuclear spins, (c) Positions 
Tfc of resonances with order k observed in (a) relative to the Larmor period Tl ~ 2n/ujL, Ar — TkjTh — (2fc — l)/4. The five 
sets of equally- spaced resonances correspond to the spins numbered in (b). Lines are fits to Eq. ^ (d) Close up for nuclear 
spin % (j ^ 8. 57. ..8. 59 //s) with A'^ = 96. Line: fit based on equation[2] Errors are ±1 standard deviation (s.d.). 
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83.8(6) 


21(1) 
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47(2) 


30(5) 


3 


55(2) 
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19(1) 


133(3) 
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33(1) 


132(1) 
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25.1(7) 


51(2) 



TABLE I: Hyperfine coupling strength LOh and angle Q for 
the six nuclear spins identified in Fig. 2. These values were 
obtained by fitting a single well-isolated resonance for each 
nuclear spin based on Eq. [2] Uncertainties are 2 s.d. 



imal coupling that can be resolved in an ESR measure- 
ment, which is given by the ESR line width. We find that 
our method detects hyperfine strengths as small as ^ 20 
kHz, about an order of magnitude smaller than the mea- 
sured line width of s/^jijiT^^) = 161(1) kHz. Further- 
more, we resolve differences in hyperfine strength down 
to - 10 kHz. 

Assuming that the interaction is purely dipole-dipole, 
the values in Table 1 correspond to distances to the NV 
center between 0.6 and 1.2 nm. The fact that we can 
distinguish multiple weakly coupled spins beyond those 
that are coupled strongest to the NV demonstrates that 
our method can be used to create tomographic images of 
the spin environment at the single nuclear spin level. 

We validate our approach by calculating the signal ex- 
pected from the values in Table 1, and comparing the re- 
sult with independent measurements over a broad range 




FIG. 3: Comparison of the measured signal with the predic- 
tion based on the parameters in Table 1 (orange line). We 
observe good agreement for the positions and amplitudes of 
multiple resonances for magnetic fields of both Ba =401 G 
and _Bo = 505 G. Error bars are ±1 s.d. 



of free evolution times at two different magnetic fields 
(Fig. 3). We find excellent agreement for both the po- 
sitions and amplitudes of the resonances, confirming the 
accuracy of the theoretical model and the determined pa- 
rameters. 

Finally, we demonstrate that we can coherently rotate 
a weakly coupled nuclear spin over a desired angle by 
tuning the number of pulses N . Figure 4(a) plots the 
signal for a selected resonance (fc = 8) of spin 3 for dif- 
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FIG. 4: Coherent conditional rotations of weakly-coupled 
nuclear spin 3. (a) P^, with diflterent numbers of pulses A'^ 
for the resonance centered at r « 8.4 /is. Lines: calculations 
based on the values in Table 1. Error bars: ±1 s.d. (b) 
at resonance as function of A^. The experimental values are 
obtained from Gaussian fits. A fit using the values in Table 1 
(line) yields a decay constant for the oscillation of 1.8(3) ms 
(supplementary information). Uncertainties are 2 s.d. 



ferent number of pulses N . The depth of the resonance 
first increases with N until the maximum contrast is ob- 
tained for TV = 7r/0 w 56. For more pulses the depth 
decreases again. In Fig. 4(b) we plot the signal at the 
center of the resonance as a function of N . 

The oscillation observed in Fig. 4(b) demonstrates the 
coherent conditional rotation of a single weakly coupled 
^■^C spin. For iV = 28 the signal reaches ~ 0.5. Here, 
the nuclear spin has rotated over an angle N(f)/2 « 7r/2, 
in a direction which is conditioned by the electron spin 
state (similar to the case illustrated in Fig. 1(c)). This 
sequence corresponds to a maximally entangling opera- 
tion, equivalent to the quantum controUed-NOT gate up 
to single-qubit rotations. For — 56, the nuclear spin 
has rotated over an angle N(j)/2 sa tt. Here, the two con- 
ditional rotations lead to the same final nuclear spin state 
up to a 27r phase difference. This phase difference trans- 
fers to the electron spin, yielding the pure state | — x) 
and signal P^, « 0. 

Unconditional rotations of the nuclear spin can be 
implemented by using different values for r (see e.g. 
Fig. 1(d)) [53]. A combination of conditional and un- 
conditional operations can be used to initialize the nu- 
clear spin by swapping its state with the electron |15j or 
for reading out the nuclear spin state in a single-shot by 
mapping it onto the electron spin [THl IHj . Our re- 
sults thus indicate the possibility of using weakly coupled 
nuclear spins as fully controllable qubits. 



The oscillation in Fig. 4(b) is damped on a timescale 
of a few ms. This timescale is consistent with the longi- 
tudinal relaxation of the electron spin at room temper- 
ature (Ti process) [55]. At cryogenic temperatures this 
relaxation time exceeds seconds [5S] , allowing for the im- 
plementation of high-precision quantum gates on weakly 
coupled -^^C nuclei. 

In conclusion, we have isolated, characterized and co- 
herently controlled individual weakly coupled nuclear 
spins embedded in a spin bath. Because we address spins 
beyond the few nearest to the NV center, our method can 
enable the tomography of ensembles of spins in diamond 
and, potentially, in external samples |29j . In addition, 
the method enables coherent gates between the electron 
spin and weakly-coupled nuclear spins. Our results thus 
indicate a clear pathway for using weakly coupled nuclear 
spins as qubits, thereby eliminating the need for strong 
coupling and greatly extending the possible number of 
qubits within a local register. 

Author's note - While finalizing this manuscript we 
became aware of two complementary studies that con- 
sider the sensing of weakly-coupled nuclear spins in the 
low magnetic field regime |30j and in isotopically purified 
diamond [31.. 
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SUPPLEMENTAL MATERIAL 



as 



SYSTEM DYNAMICS DURING THE 
DECOUPLING SEQUENCE 

In this section we derive the equations used in the main 
text. We consider the conditional dynamics of a nuclear 
spin that interacts with the NV electron spin through the 
hyperfine coupling. 



Hamiltonian of the system 

With an appropriate rotation of the coordinate axes, 
the Hamiltonian of the NV spin coupled to a single ^'^C 
spin is given by 

H = ASJ, + BSJI+ujlIz = |0)(0|i/o + |l)(l|ffi, (4) 

where Si (li) are the Pauli matrices of the electron (nu- 
clear) spin, ujl is the nuclear Larmor frequency, and 
A = uJhCOsO [B = UhSinO) is the parallel (perpen- 
dicular) component of the hyperfine coupling (see Fig. 
1(a)). The nuclear spin evolves according to a Hamilto- 
nian which is conditioned by the electron spin state: Hq 
if the electron is in rag — (state |0)), and Hi if the 
electron is in = — 1 (state |1)), with 

Ho ^ lolIz, and Hi = {A + ujl)Iz + BI^. (5) 



Net result of the decoupling sequence 

As described in the main text, we probe the conditional 
interaction by preparing the electron spin in \x) — (|0) -I- 
|l))/-\/2, driving it by a decoupling sequence with N/2 
decoupling units of the form t — tt — 2t — tt — r (Fig. 1 (b) ) , 
and measuring the x-component of the NV electron spin 
at the end of the sequence. 

The probability that the initial state \x) is preserved 
is given by = (M + l)/2, with 



M = Re Tr(;7o?7i ) 



(6) 



where Uq and Ui are the conditional evolution operators 
for the nuclear spin corresponding to the initial states of 
the NV spin = and TOs = — 1 respectively. These 
operators, describing the evolution after N/2 decoupling 

- N/2 - N/2 

units, have the form Uq — Vq and Ui — Vi , where 
Vb and Vi are the conditional evolution operators for a 
single decoupling unit {N = 2): 

Vq = exTp[—iHoT]exp[—iHi2T]exp[—iHQT] (7) 
Vi — exp[—iHiT]exp[~iHo2T]exp[—iHiT]. (8) 

Since any unitary evolution of a single qubit is a rotation, 
the conditional operators Vq and Vi can be represented 



expl 
exp[ 



(I • fio) 
(I-fii) 



(9) 
(10) 



which illustrates that the net evolution of the nuclear spin 
after a single decoupling unit is a rotation by an angle 
(j) around an axis that depends on the initial state of 
the electron spin: no for initial state = and fii for 
initial state mj — —1. Note that the rotation angle 
is independent of the electron spin input state because 
COS0 = TvVn = TvVi. 

With these expressions, we obtain Eq. 2 of the main 
text: 



M = 1 - (1 - fio • rii) sin^ 



N(j) 



(11) 



Thus, we only need to determine the net rotation angle (p 
and the inner product no • rii to completely characterize 
the signal. 



Nuclear spin rotation axes and angle 



Using eqs.[7]-[T0j we obtain expressions for the rotation 
angle and the inner product fio • fii : 



1 - fio • fi 



COS0 = cos a cos /3 — TOz sin a sin /? (12) 
2 (1 — cosa)(l — cos/3) 



1 4- cos a cos j3 — ruz sin a sin 



5(13) 



where — {A -\- ljl)/Cj, = B/uj, and to — 

{A + lol)^ + B^. For simplicity, we introduced the an- 
gles a = UT and (3 — lolT- 



High magnetic field: 2> ijJh 

The advantage of performing the experiments at high 
magnetic field, so that oj^ ^ w/j, is that the signal M has 
a clear form: for most values of r, the depth of the signal 
modulation (1 — fio-fii) is small, of order of to^ ^ B^ /uj\. 
Here, the axes fio and fii are practically parallel and the 
nuclear spin undergoes an unconditional rotation (Fig 
1(d)). The modulation is only strong when the decou- 
pling sequence is precisely resonant with the nuclear spin 
dynamics. At resonance, the axes fio and fii are anti- 
parallel, i.e. fio • fii = —1, leading to a conditional ro- 
tation of the nuclear spin (Fig. 1(c)). The condition for 
the resonances is 



a 13 I 
tan — tan — = — 
2 2 rriz 



(14) 



For the high magnetic field case where <C 1, in 
first order in this small quantity, the solution for Eq. [14] 
is a = — (/? -|- tt) -|- 2fc7r with any integer fc, so that 
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tan (a/2) = l/tan(/?/2), and the corrections are of the 
second order in rux- In this case the inter-pulse distance 
2t corresponding to the resonance condition satisfies 



(2/c + l)7r 
2wL + A ' 



(15) 



which is equivalent to the Eq. 3 of the main text. 

We now provide more insight into the shape and am- 
plitude of the resonances. To understand what happens 
when r is varied from the resonant value r^, we assume 
that T = Tfe(l-I-A) with small A. Then, the inner product 
of the rotation axes becomes 



l-{no-n^)^2/{l + 5l/ml), 



(16) 



where 5^ = {2k + l)7rA, and we assume as above that 
v4, i? <C uji^. Thus, this modulation depth as a function 
of A is a Lorentzian with a small width of mx/[{2k + 
l)7r] (this justifies our assumption about A being small, 
not exceeding rrix ^ 1 too much). The final result as a 
function of 5k acquires the form 



M = 1 



1 



5llml 



N 



SI 



(17) 



This expression demonstrates that even when mx <^ 1 
the resonance can be visible. When the argument of the 
sine function above is small, the amplitude of the res- 
onance is of order m^N'^, and is therefore amplified by 
a factor of N'^ . By increasing the number of pulses, we 
increase the sensitivity of the detection of the nuclear 
spins, i.e. we can address more and more nuclei, which 
are coupled weaker and weaker to the NV spin. 



ELECTRON SPIN RELAXATION AND 
DEPHASING 



For the fits in the main text we take the overall de- 
cay of the signal due to the relaxation and dephasing 
of the electron spin into account. The equation relating 
the obtained signal to the modulation M due to the 
interaction with the nuclear spins becomes: 



P. 



l/2e-2JVr/T„^^ + 1/3 + l/6e-2^^/^^, 



(18) 



in which Ta and Tp are two phenomenological time con- 
stants. Note that at r = the signal modulation is 
maximum and equation [18] is equivalent to the ideal case 
given in the text. 

describes the Ti relaxation of the electron spin into 
the rris = +1 sublevel. This relaxation affects the center 
value that the modulation M acts on. We measure this 
time constant by taking the average value of the outcomes 
of measurements along X and —a; (| — a;) = (|0) — |1))/a/2), 
i.e. effectively we measure for M — 0. We obtain = 
2.2 ms. 

accounts for dephasing of the electron spin and in- 
cludes not-decoupled dephasing, microwave pulse errors 
and decoherence due to Ti relaxation. These effects re- 
duce the total modulation amplitude. This reduction is 
only approximately exponential, but this is a good ap- 
proximation over the small r ranges that are used for 
the fits. The values obtained from fits vary between 1.3 
and 2.2 ms for different measurements depending on the 
number of pulses, the fidelity of the 7r-pulses, and the 
applied magnetic field. 



